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Abstract 

We report a numerical investigation of surface plasmon (SP) propagation in ordered and dis- 
ordered linear chains of metal nanospheres. In our simulations, SPs are excited at one end of a 
chain by a near-field tip. We then find numerically the SP amplitude as a function of propagation 
distance. Two types of SPs are discovered. The first SP, which we call the ordinary or quasistatic, 
is mediated by short-range, near-field electromagnetic interaction in the chain. This excitation is 
strongly affected by Ohmic losses in the metal and by disorder in the chain. These two effects re- 
sult in spatial decay of the quasistatic SP by means of absorptive and radiative losses, respectively. 
The second SP is mediated by longer range, far- field interaction of nanospheres. We refer to this 
SP as the extraordinary or non-quasistatic. The non-quasistatic SP can not be effectively excited 
by a near-field probe due to the small integral weight of the associated spectral line. Because of 
that, at small propagation distances, this SP is dominated by the quasistatic SP. However, the 
non-quasistatic SP is affected by Ohmic and radiative losses to a much smaller extent than the 
quasistatic one. Because of that, the non-quasistatic SP becomes dominant sufficiently far from 
the exciting tip and can propagate with little further losses of energy to remarkable distances. 
The unique physical properties of the non-quasistatic SP can be utilized in all-optical integrated 
photonic systems. 
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I. INTRODUCTION 



Surface plasmons (SPs) are states of polarization that can propagate along metal- 
dielectric interfaces or along other structures without radiative losses. Polarization in an 
SP excitation can be spatially confined on scales that are much smaller than the free-space 
wavelength. This property proved to be extremely valuable for manipulation of light en- 
ergy on subwavelength scales [1, 2], miniaturization of optical elements [3], and achieving 
coherent temporal control at remarkably short times [4, 5]. SP excitations in ordered one- 
dimensional arrays of nanoparticles have attracted significant attention in recent years due 
to numerous potential application in nanoplasmonics [6-10]. A periodic chain of high con- 
ductivity metal nanospheres can be used as an SP wave guide - an analog of an optical 
waveguide [11]. High-quality SP modes in ordered and disordered chains may be utilized 
in random lasers [8]. Electromagnetic forces acting on linear chains of nanoparticles can 
produce the effect of optical trapping [12]. Various spectroscopic and sensing applications 
have also been discussed [13-15]. 

In this paper we study theoretically and numerically propagation of SP excitations in 
long ordered and disordered chains of nanospheres. Although, under ideal conditions, SP 
excitations can propagate without loss of energy, in practical situations this is not so. There 
are two physical effects that can result in decay of SP excitations as they propagate along 
the chain. The first effect is Ohmic losses due to the finite conductivity of the metal. The 
second effect is radiative losses due to disorder in the chain (scattering from imperfections) . 
This effect is more subtle and is closely related to the phenomenon of locahzation. In this 
paper we discuss both effects and illustrate them with numerical examples. 

We first focus on decay due to Ohmic losses and show that it can be suppressed at 
sufficiently large propagation distances. The main idea is based on exploiting an exotic 
non-Lorentzian resonance in the chain which originates due to radiation-zone interaction 
of nanoparticles [13, 16] and can not be understood within the quasistatics, even when 
both the nanoparticles in the chain and the inter-particle spacing are much smaller than 
the wavelength. Prom the spectroscopic point of view, the non-Lorentzian resonances are 
manifested by very narrow lines in extinction spectra [13, 14]. One of the authors (V.A.M.) 
has argued previously that the small integral weight of the spectral lines associated with these 
resonances precludes them from being excited by a near- field probe [16] . This property would 
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make the non-Lorentzian resonance a curiosity which is rather useless for nanoplasmonics. 
However, numerical simulations shown below reveal that the corresponding SP has relatively 
small yet nonzero amplitude and is also characterized by very slow spatial decay. Therefore, 
in sufficiently long chains, this SP becomes dominant and can propagate, without significant 
further losses, to remarkable distance. Wc stress that the non-Lorentzian SP is an excitation 
specific to discrete systems; it does not exist, for example, in metal nanowires. 

But the above consideration applies only to ordered chains. Therefore, we consider next 
the effects of disorder. To isolate radiative losses due to scattering on imperfections from 
Ohmic losses, we consider nanoparticles with infinite conductivity (equivalently, zero Drude 
relaxation constant). Although such metals do not exist in nature, an equivalent system 
can be constructed experimentally by embedding metallic particles into a dielectric medium 
with positive gain. We show that while the ordinary (defined more precisely in the text 
below) SP excitations are very sensitive to off-diagonal (position) disorder, the SP due to 
the non-Lorentzian resonance is not. Diagonal disorder (disorder in nanoparticle properties) 
is also considered. 

The paper is organized as follows. In section II, we describe the theoretical model and 
introduce basic equations. Conditions for resonance excitation of SPs in a chain are con- 
sidered ins section III. Numerical results for propagation in ordered and disordered chains 
are reported in sections IV and V, respectively. Finally, section VI contains a summary of 
obtained results. 

II. THEORETICAL MODEL 

Consider a linear chain of nanosphercs with radiuses a„ centered at points x„. We 
work in the dipole approximation which is valid if Xn+i — (on+i + an)/2 and has been 
widely used in the literature [7, 8, 10, 17-19]. The n-th nanosphere is then characterized by a 
dipole moment with amplitude dn oscillating at the electromagnetic frequency cu. The dipole 
moments are coupled to each other and to external field by the coupled-dipole equation [13] 



where q;„ is the polarizability of the n-th nanosphere. En is the external electric field at 
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the point x„, k — uj/c is the free space wave number and Gk{x,x') is the appropriate 
element of the free space, frequency-domain Green's tensor for the electric field. The latter is 
translationally invariant with respect to spatial variables, namely, Gk{x, x') = Gk{x—x', 0) = 
Gk{x' — a;,0). For an SP polarized perpendicular and parallel to the chain, the respective 
functions G^ and G\ are given by 
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where a^n^^ is the Lorenz-Lorentz quasistatic polarizability of a sphere of radius a„ and 
2ik^/3 is the first non-vanishing radiative correction to the inverse polarizability; account of 
this correction is important to ensure that the system conserves energy [20]. The Lorenz- 
Lorentz polarizability is given, in terms of the complex permeability of the n-th nanosphere 
en, by 
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We further adopt, for simplicity, the Drude model for e„: 
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where uj is the electromagnetic frequency, uj^n is the plasma frequency, and 7„ is the Drude 
relaxation constant in the n-th nanosphere. The inverse polarizability of the n-th nanosphere 
is then given by 
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/-s/3 is the Prohlich frequency. The SP resonance of an isolated n-th 



nanosphere takes place when a; = a;Fn- Polarizability q;„ at the Prohlich resonance is purely 
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imaginary; if, in addition, there are no Ohmic losses in the material (7^ = 0), the resonance 
value of the polarizability becomes = ctres = —iS/2k^, irrespectively of the particle radius. 

Suppose that SP is excited at a given site (say, n = Uq) by a near- field probe. Then the 
external field can be set to En = -Eo^n.no- C)f course, this is an idealization: the field produced 
even by a very small near-field tip is, strictly speaking, non-zero at all sites. However, this 
approximation is physically reasonable because of the fast (cubic) spatial decay of the dipole 
field in the near-field zone. The solution with En — 5n,no is, essentially, the Green's function 
for polarization. We denote this Green's function by 'Dfc(x„, It satisfies 

where either (2) or (3) should be used for G^, depending on polarization of the SP. In the 
case of a finite or disordered chain, one can find 'Dfc(,x„, by solving Eq. (1) numerically. 
However, in infinite ordered chains such that a„ = a = const and Xn+i — Xn = h = const, 
the following analytic solution is obtained by Fourier transform [16]: 

-n n\ exp(ig3;„) hdq 

J_^/h 1/a - S{k,q) 2tt 

where S{k, q) is the "dipole sum" given by 

S(k, q)^2j2 Gk{0, Xn) cos(gx„) . (10) 

n>0 

Obviously, in infinite chains ^^^{xn, Xn') — T^ki^n ~ ^n', 0) = ^^^{xn' — 0). Note that the 
dipole sum (10) is independent of material properties. It can be shown [13] that, for all values 
of parameters, Im5'(A;, q) > —2k^/3. The equality holds when q > k. This is a manifestation 
of the fact that SPs with q > k are non-radiating due to the light-cone constraint [8]. 
Non-radiating modes exist if n/h > k or, equivalently, if A > 2h. Obviously, these SPs do 
not couple to running waves but can be excited by a near-field probe. The dimensionless 
radiative relaxation parameter can be defined as Q{k, q) — \hn.S{k, q) -\-2k^ /^]/{2k^ /?,)■, this 
factor is identically zero ior q> k. 
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III. DISPERSION RELATIONS AND RESONANT EXCITATION OF SP 



In this section, we consider periodic chains with an — a, Upn — Up, — ojf, — 1 
and, consequently, with the constant polarizabihty = a. 

It follows from formula (9) that the wave numbers q of SP excitations that can propagate 
effectively in an infinite periodic chain are such that l/a ~ S{k,q). This condition must 
be satisfied for some range of q which is, in some sense, small. Indeed, there is no effective 
interaction in the chain if l/a — S{k, q) ~ const, in which case integration according to (9) 
yields [21] (i„ oc 5nfi- We now examine the conditions under which the denominator of (9) can 
become small. The dispersion relation in the usual sense, e.g., the dependence of the resonant 
SP frequency on its wavenumber is obtained by solving the equation l/a{uj) —S{uj/c, q) — 0. 
This approach was adopted, for example, in Refs. 7, 8, 10, 17. Solution to the above 
equation depends on the model for e{uj) and may result in several branches of the complex 
function u!{q). Here we adopt a slightly different point of view. Namely, we note that for 
real frequencies cu and wave numbers q, the real part of the denominator can change sign 
while the imaginary part is always non- negative. Physically, Re[l/Q;(a;) — S{u!/c,q)] can 
be interpreted as the generalized detuning from a resonance while Im[l/Q;(a;) — S{u!/c,q)] 
gives total (radiative and absorptive) losses. We thus define the resonance condition to be 
Re[l/a — S{k, q)] = and view uj and q as independent purely real variables. 

Plots of the dimensionless function h^S{k, q) are shown in Fig. 1 for some typical sets of 
parameters and for two orthogonal polarizations of the SP. Note that the imaginary part 
of h^S{k,q) is related to radiative relaxation parameter Q{k,q) [shown in Fig. l(c,d)] by 
hHm.S{k,q) — [2(A;/i)^/3][(5(A;, g) — 1]. Apart from the very narrow peaks appearing in 
the case of perpendicular polarization and centered at g = /c, [22] the numerical values of 
/i^lReS*! do not exceed, approximately, 6. On the other hand, according to (7), we have 
h?'Re{l/a) = {h/aY[l — (cu/o^p)^]. The dipole approximation is vafid when {h/aY > 64. 
Therefore, if we stay within the range of parameters in which the dipole approximation is 
valid, the real part of the denominator in (9) is relatively small only ii cu ~ cup, i.e., if cu is 
near the Frohlich frequency of an isolated sphere. This is the case that will be considered 
below. 

Consider first oscillations polarized orthogonally to the chain and let uj = ujp. The 
condition of resonant excitation of SP is then ReS{k, q) = 0. It can be seen from Fig. 1(a) 
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that, for sufficiently small values of the dimensionless parameter kh^ there are two different 
values of q that satisfy the above resonance condition. 

The first solution is g = gi « 0.477r//i (for kh = 0.27r). The value of gi depends only 
weakly on kh^ as long as kh < 0.27r, and is approximately the same as in the quasistatic 
limit kh = 0, in which case qi ^ 0.467r//i. We will refer to the SP with this wave number as 
the ordinary, or the quasistatic SP. This is because propagation of this SP is mediated by 
near-field interaction, while the far-field interaction is suppressed by destructive interference 
of waves radiated by different nanoparticles in the chain. Since the wave number of the 
ordinary SP is greater than k, it propagates without radiative losses. However, it may 
experience spatial decay due to absorption in metal. The characteristic exponential length 
of decay can be easily inferred from (9) by making the quasi-particle pole approximation. 
Namely, we approximate S{k, q) as 



S{k,q) p^ReS{k,qi) + {q - qi) 



dRfiS{k, q) 



2k 



3 



-i^. (11) 
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dq 

We then extend integration in (9) to the real axis and obtain the following characteristic 
exponential decay length i: 



where 



dReSjk, q) 
dq 



(12) 

9=91 



5 = -Im(l/Q;) -2A;V3 (13) 

is a positive parameter characterizing the absorption strength of the nanosphere. In general, 
it can be shown that (5 = in non-absorbing particles whose dielectric function eiuj) is 
purely real at the given frequency uj. For the Drude model adopted in this paper, we have 
5 — uj^/a^uj'^. Thus, the ordinary SP can decay exponentially due to absorption in metal 
with the characteristic scale given by (12). 

The second solution is obtained at g = g2 ~ ^, when ReS'(A:, q) has a narrow sharp peak 
as a function of q (for fixed k) [23]. This peak is explained by far-field interaction in an 
infinite chain [13, 16]. It does not disappear in the limit kh 0, but becomes increasingly 
(super-exponentially) narrow [16]. Note that in the above limit, this peak appears as a 
singularity of zero integral weight which can not be obtained from the formal quasistatic 



7 




FIG. 1: (a,b): Real parts of the dipole sum S{k,q) as a function of q, for different values of 
kh, as indicated. Sharp peaks corresponding to divergence of JieS{k,q) at the point q = k are 
not completely resolved. Quasistatic results obtained in the limit k = are shown by curves 
labeled "QS". Only half of the Brillouin zone is shown since S{k,—q) = S{k,q). Polarization of 
SP is perpendicular (a) and parallel (b) to the chain. (c,d): Dimensionless radiative relaxation 
parameter Q{k,q) = [lmS{k,q) + 2k^ /3]/ {2k^ /3) function of q for the same sets of parameters 
as above. Polarization of SP is perpendicular (c) and parallel (d) to the chain. Quasistatic result 
is not shown since Q is not defined in the quasistatic limit. 
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approximation. We will refer to this SP as extraordinary, or non-quasistatic. It is mediated 
by the far-field interaction. The latter is important in the case of extraordinary SP because 
of the constructive interference of far field contributions from all nanoparticles arriving at 
a given one. Obviously, the extraordinary SP can be excited only in chains with h < A/2, 
where A is the wavelength of light in free space. 

It is interesting to consider radiative losses of the extraordinary SP. As can be seen from 
Fig. 1(c), the factor Q{k,q) is discontinuous at q — k: it is zero ior q > k but positive 
for q < k. Since the extraordinary SP has q ~ k, it can experience some radiative losses, 
although the exact law of its decay depends on parameters of the problem in a complicated 
manner. This is, in part, related to inapplicability of the quasi-particle pole approximation 
for evaluating the integral (9) in the vicinity of g ~ g2- 

In the case of oscillations polarized along the chain, the resonance condition can be 
satisfied only at q — qi ^ OABir/h. This is the wave number of an ordinary (quasistatic) 
SP which depends on k only weakly, as long as kh < 0.27r. However, the extraordinary 
(non-quasistatic) SP can be excited even for longitudinal oscillations. Mathematically, this 
can be explained by observing that dReS{k,q)/dq diverges at q = k while dlmS{k,q)/dq 
is discontinuous at q = k and performing integration (9) by parts. However, the amplitude 
of the extraordinary SP is much smaller for longitudinal oscillations than for transverse 
oscillations; this will be illustrated numerically in the next section. 

Finally, we note, that when inter-sphere separations become smaller than the radiuses, 
higher-order multipole resonances can be excited [24-26]. In this case, resonant excitation 
of SP can become possible even at frequencies which are far from the Frohhch resonance of 
an isolated sphere, e.g., in the IR part of the spectrum. 

IV. PROPAGATION IN FINITE ORDERED CHAINS 

We now turn to propagation of SP in ordered chains of finite length N. We, however, 
emphasize that the finite size effects play a very minor role in the computations shown below. 
Citrin [18] has studied dispersion relations in finite chains and has found that the infinite- 
chain limit is reached at ~ 10 (although wc anticipate that longer chains arc needed for 
accurate description of the extraordinary SP). In this and following sections, we work with 
chains oi N > 1000. In this limit, propagation of both ordinary and extraordinary SP is not 
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much different from the case of infinite chains. In particular, we have verified numerically 
that the Green's function X'fe(x„, X501) (for kh — 0.27r) in a chain oi N — 1001 particles 
does not differ in any significant way from that in an infinite chain, except for values of n 
very close to either end of the finite chain. The Green's function Gk{xn,xi) (here Uq = 1 
is the end-point of the finite chain) differed by a trivial factor in finite and infinite chains 
(results not shown). However, proper numerical evaluation of integral (9) required very 
fine discretization of q and was a more demanding and less stable procedure than direct 
numerical solution of the system of equations (8). 

In this section, we take Un — a — const, an — en — const and x„ — nh, n — 1, . . . ,N. 
We also assume that kh = 0.27r and h = 4a. Practically, this can be realized for silver 
particles in a transparent host matrix with refractive index of approximately n = 1.4 so that 
the wavelength at the Frohlich frequency is Ap — 2-Kc/uj-p ~ 400nm, the chain spacing is 
h — 40nm = O.lAp and the sphere radius is a = lOnm = /i/4. The dipole approximation is 
very accurate for this set of parameters. We then obtain the SP Green's function T>k{xm xi) 
by solving Eq. (8) numerically. 

The absolute value of the normalized SP Green's function 

in a chain oi N — 1000 nanospheres is shown in Fig. 2 as a function of x (sampled at x — x„) 
for two orthogonal polarizations. Here the frequency of SP was taken to be exactly equal 
to the Frohlich frequency and the Drudc relaxation constant was 7 = 0.002a;F- It can 
be seen from the figure that two different SP are excited in the system. The first is the 
ordinary (quasistatic) SP that decays exponentially as exp{—x/£), where i is defined by 
(12). Corresponding asymptotes are shown by dotted fines. Note that the quasistatic SP 
in a finite chain (with the point of excitation coinciding with one of the chain ends) is very 
well described by the exponential decay formula, even though the latter was obtained for 
infinite chains. When the amplitude of the ordinary SP becomes sufficiently small, there is 
a crossover to the extraordinary SP. The decay rate of the extraordinary SP is much slower. 
We have also confirmed by inspecting the real and imaginary parts of 'Dk{xn,Xi) (data not 
shown) that it oscillates at the spatial frequency corresponding to the ordinary SP in the 
fast-decaying segments of the curves shown in Fig. 2 and with the spatial frequency that 
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FIG. 2: Propagation of a SP in an ordered chain of = 1000 nanospheres for orthogonal (ORT) 
and parahel (PAR) polarization of oscillations with respect to the chain. Parameters: uj = lof, 
7/tJF = 0.002, A = Wh, h = 4a. 

corresponds to the extraordinary SP in the slow-decaying segments. 

Mathematically, the relatively slow decay of the extraordinary SP can be understood as 
follows. First, note the the exponential decay of the ordinary SP is, in fact, the result of 
superposition of an infinite number of plane waves whose wave numbers are in the interval 
Aq oc 6. The corresponding wave packet decays spatially on scales i oc 1/Ag. However, in 
the case of extraordinary SP, Aq can not be defined since the corresponding resonance is 
non-Lorentizian. It can be, however, stated that the extraordinary SP is a superposition of 
plane waves whose wave numbers are very close to k. Since, however, the wavenumbers can 
still slightly deviate from k, some spatial decay at large distances can still occur. 

The conclusion we can make so far is that the ordinary SP experiences exponential decay 
along the chain due to Ohmic losses in metal. This decay is very accurately described by 
the quasi-particle pole approximation. The extraordinary SP has, initially, much smaller 
amplitude than the ordinary one. This is because the peaks in Fig. 1 are very narrow. The 
quasi-particle pole approximation is invalid for the extraordinary SP and its decay is much 
less affected by Ohmic losses. As a result, the extraordinary SP decays at a much slower rate 
and, at sufficiently large propagation distance, begins to dominate. We also note that the 
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FIG. 3: Same as in Fig. 3 for different ratios "y/iv-p and for SP polarized orthogonally (a) and 
parallel (b) to the chain. 

extraordinary SP can be excited even for longitudinally-polarized SP, although its amplitude 
is smaller by some four orders of magnitude than for the case of transverse oscillations. 

In Fig. 3, we illustrate the influence of Ohmic losses on SP propagation. Here we plot 
as a function of x (sampled at x = Xn) for uo = uj-p and different values of the ratio 
/uj-p. First, in the absence of absorption (7 = 0), the ordinary SP propagates along the chain 
without decay. Once we introduce absorption, the ordinary SP decays exponentially with 
the characteristic length scale i given by (12). Note that, for the specific metal permeability 
model (6), £ oc ujp/'j. Some dependence of the rate of decay of the extraordinary SP on 
the ratio 'y/oop is visible in the case of orthogonal polarization [Fig. 3(a)]. However, when 
the polarization is longitudinal [Fig. 3(b)], decay of the extraordinary SP is dominated by 
radiative losses. In particular, the slow-decaying segments of the curves for j/up = 0.002 
and j/ujp = 0.004 in Fig. 3(b) coincide with high precision. 

Next, we study SP propagation for different values of the ratio uj/ujp. As noted above, 
we assume the parameters kh = uoh/c = 0.27r and h/a = 4 to he fixed. Thus, u/up can 
vary either due to a change in ujp or due to a simultaneous change in u, h and a such 
that uh = const and h/a = const. It follows from Fig. 1(a) that, for the selected set of 
parameters, the ordinary plasmon can be excited for —0.89 < /i^Re(l/a) < 1.56. This 



12 




FIG. 4: Same as in Fig. 3 for different ratios uj/to-p and 7/wF) as indicated, and for SP polarized 
orthogonally to the chain. 

corresponds to oj/ojy lying in the interval 0.988 < cj/c^f < 1.007. In Fig. 4, we illustrate 
propagation of SP excitations for some values of uj/uj-p inside this interval, exactly at the 
lower and upper bounds of this interval, and slightly outside of the interval. Results are 
shown for two values of absorption strength: 7/tt;F = and 7/cc'F = 0.002. 

First, consider the two cases when uj/uo-p is outside of the interval where the ordinary 
SP can be excited: oj/up = 0.984 and u/up = 1.010 (correspondingly, h^Re{l/a) = 2.03 
and /i^Re(l/a) = —1.29), shown in Fig. 4(a,f). In the case u/up = 0.984, the ordinary 
SP exhibits very fast spatial decay, which is characteristic for the noninteracting limit when 
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Gk{n, no) oc 5n,no- After the initial decay of the ordinary SP, the extraordinary SP becomes 
dominating. The extraordinary SP decays slowly by means of radiative losses. It is inter- 
esting to note that decay of the extraordinary SP is almost unaffected by Ohmic losses in 
metal, although a noticeable dependence on 7 appears if we further increase the parameter 
7/cjF by the factor of 10 (data not shown). A qualitatively similar behavior is obtained 
at lu/luf — 1.010. However, the decay of extraordinary SP in this case is a little faster. 
Paradoxically, the curve corresponding to ^/lof — 0.002 is slightly higher than the curve 
corresponding to 7 = in Fig. 4(f) (similar pecuharity is seen in Fig. 4(e)). 

When the ratio cu/cuf is inside the interval where the ordinary SP can be excited 
[Fig. 4(c,d)], SP propagation is strongly influenced by absorptive losses. In the absence 
of such losses, the ordinary SP propagates along the chain indefinitely and dominates the 
extraordinary SP. However, in the presence of even small absorption, the ordinary SP de- 
cays exponentially so that, at sufficiently large propagation distances, the extraordinary SP 
starts to dominate. A qualitatively similar picture is also obtained for the borderline case 
u/uf — 0.984 [Fig. 4(b)]. In the second borderhne case [Fig. 4(e)], SP propagation is more 
complicated. The wave numbers of both ordinary and extraordinary SPs in this case are 
close to k, so that both can experience radiative decay, as is evident in the case of zero 
absorption. 

To conclude this section, we note that propagation of SP excitations in long periodic 
chains can be characterized by exponential decay. This decay is caused either by absorptive 
or by radiative losses. At small propagation distances, energy is transported by the ordinary 
SP excitation, if the ordinary SP can be excited (e.g., if Re(l/Q;) is inside the appropriate 
interval). However, at sufficiently large propagation distances, there is a cross over to trans- 
port by means of the extraordinary SP. In this case, propagation is mediated by far-zone 
interaction and is characterized by slow, radiative decay which is affected by absorptive 
losses only weakly. 

We note that exponential decay in ordered chains, if exists, is not caused by Anderson 
localization, since we have not, so far, introduced disorder into the system. For example, as 
was discussed in section 111, a linear superposition of delocalized plane wave modes of the 
form (9) can exhibit exponential decay with the characteristic length (12). The irreversible 
exponential decay is, in fact, obtained because the delocalized modes form a truly continuous 
spectrum (are indexed by a continuous variable q) . Of course, any superposition of discrete 
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delocalized modes would result in Poincare recurrences. 

V. PROPAGATION IN DISORDERED CHAINS 

The ordinary (quasistatic) SP propagates in ordered chains without radiative losses due to 
the perfect periodicity of the lattice. However, once this periodicity is broken, the quasistatic 
SP can experience radiative losses and spatial decay even in the absence of absorption. 
Dependence of the radiative quality of SP modes in finite one-dimensional chains on disorder 
strength was studied in Ref. 8. It was shown that position disorder tends to decrease the 
radiative quality factor of initially non-radiating ( "bound" ) modes. In this section we study 
how the disorder influences propagation of the SP along the chain and take a separate look 
at the ordinary and extraordinary SPs. We also consider two types of disorder: off-diagonal 
and diagonal. Off-diagonal disorder is disorder in particle positions while all particles are 
identical. Diagonal disorder arises due to differences in particle properties, even if the particle 
positions are perfectly ordered. 

The exponential decay due to Ohmic losses in the material can mask the effects of disorder. 
Therefore, we assume in this section that the nanospheres are non-absorbing, i.e., set 7^ — 0. 
Physically, absence of absorption can be realized by embedding the chain of nanoparticles in 
a transparent dielectric host medium with positive gain [27-30]. Such medium has dielectric 
permeability eh with positive real and negative imaginary parts. The gain can be tuned so 
that the effective permeability of inclusions e/eh is purely real and negative. We also work 
in the regime kh — 0.27r, h — 4a. 

A. Off-Diagonal Disorder 

Off-diagonal disorder is disorder in particle position. It is called "off-diagonal" because it 
affects only off-diagonal elements of the interaction matrix In the simulations 

shown below, coordinates of particles in a disordered chain were taken to he Xn — h{n + 
where ^„ is a random variable evenly distributed in the interval [—A, A]. The random num- 
bers a-re taken to be mathematically independent. Therefore, the disorder is uncorrelated. 

The first example concerns propagation in an non-absorbing chain of N = 10, 000 
nanospheres excited exactly at the Frohlich frequency uu = uup. Numerical results for dif- 
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ferent levels of disorder are illustrated in Fig. 5 where we plot as a function of x 

(sampled aX, x — Xn) ■ 

One obvious conclusion that can be made from inspection of Fig. 5 is that disorder causes 
spatial decay of SP. Since the system has no absorption, energy is lost to radiation. However, 
the exact law of decay strongly depends on particular realization of disorder. In Figs. 6 and 7, 
we plot the function for A — 0.01 and A — 0.02, respectively, and for three different 

realizations of disorder (without the use of logarithmic scale). Giant fluctuations in the 
amplitude of transmitted SP are quite apparent. In all cases, after some initial growth, the 
amplitude decays, although some random recurrences (due to re-excitation) can take place. 
However, the amplitude |jFfc(a;)| at a given site x = Xn strongly depends on realization of 
disorder and can be very far from its ensemble average, even at very large propagation 
distances. Therefore, the function J-k{x) appears to be not self-averaging. In particular, any 
particular realization of the intensity 2a.(x) = |^fe(x)p does not satisfy either the radiative 
transport equation or the diffusion equation. The absence of transport and self- averaging is 
characteristic for Anderson locahzation. We note that the localization discussed here is an 
interference (radiative) effect which is different from localization of quasistatic polarization 
modes studied in Refs. 31, 32. Localization properties of the eigenmodes in disordered chains 
is a subject of separate investigation and will be reported elsewhere. 

The situation is comphcated by the presence of two types of SP excitations. One can 
argue that locahzation properties of these two types of SPs might be different. To investigate 
this possibility, we perform two types of numerical experiments. First, we repeat simulations 
illustrated in Fig. 5 but for lu/luf = 0.948. In this case, the ordinary SP is not excited in 
ordered chains (sec Fig. 4(a)). Results are shown in Fig. 8 . It can be concluded from the 
figure that off-diagonal disorder does not result in additional decay of the extraordinary SP. 
In fact, the curves with A — 0, A — 0.1 and A — 0.04 are indistinguishable. This should 
be contrasted with the case lo/lof — 1, when, at the level of disorder A — 0.02, spatial 
decay is already well manifested. The conclusion one can make is that the extraordinary 
SP excitations do not experience localization due to uncorrelated off-diagonal disorder in 
the chain. Therefore, spatial decay seen in Figs. 5 and 7 is decay of the ordinary SP. At 
very large propagation distances, J-'k{x) continues to oscillate around the baseline of the 
extraordinary SP. 

The second numerical experiment that can elucidate the influence of disorder on ordinary 
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FIG. 5: Propagation of SP in a chain of = 10, 000 non-absorbing (7 = 0) nanospheres for 
different levels of off-diagonal disorder A, as indicated. SP polarization is orthogonal to the chain. 
Other parameters: kh = 0.2it, uj/lof = 1, h/a = 4. 
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FIG. 6: Propagation of SPs in a chain with off-diagonal disorder of amplitude A = 0.01 for different 
random realizations of disorder. Other parameters same as in Fig. 5. 
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FIG. 7: Same as in Fig. 6, but for A = 0.02. 
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FIG. 8: Same as in Fig. 5 but for uj/oj-p = 0.984 and for different levels A of off-diagonal disorder, 
as indicated. Note that the curves for A = and A = 0.01 and A = 0.04 are indistinguishable. 

and extraordinary SPs is simulation of an extinction measurement when the chain is excited, 
instead of a nera-field tip, by a plane wave i?o exp(iga;). Namely, we will look at the depen- 
dence of the specific (per one nanosphere) extinction cross section on q. There are two 
different experimental setups that can be used to measure cTe(g). When q is in the interval 
< g < fc, this experiment can be carried out simply by varying the angle between the 
incident beam and the chain. However, values q > k are not accessible in this experiment. 
In this case, the chain can be placed on a dielectric substrate and excited by evanescent 
wave originating due to the total internal reflection of the incident beam. The maximum 
longitudinal wave number of the evanescent wave is nk, where n > 1 is the refractive index 
of the substrate. We note that it is not realistically possible to access all wave numbers up 
to q = it/ h in this way because, in the particular case kh = 0.27r, this would require the 
refractive index n = 5. Refractive indices of such magnitude are not achievable in the optical 
range. However, in a numerical simulation, we can assume that a hypothetical transparent 
substrate with n = 5 exists. Besides, all wave numbers q in the first Brillouin zone of the 
lattice can be accessible for a different choice of parameters (particularly, for larger values 
of kh). 

The specific extinction ae{q) is given by the following formula: 
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FIG. 9: Specific extinction cje as a function of lateral wave number of incident wave, q, for different 
levels of off-diagonal disorder, A. 



where dn is the solution to Eq. (1) with the right-hand side En = -Eq exp(igx„). In Fig. 9, 
we plot the dimensionless quantity h^'^a^.^q) as a function of q for transverse oscillations in 
partially disordered non-absorbing chains. Two peaks are clearly visible in the spectrum. 
The first peak at g = gi ~ OATir/h corresponds to excitation of the ordinary SP. The second 
peak at g = g2 ~ ^ corresponds to excitation of the extraordinary SP. In infinite, periodic 
and non-absorbing chains, the spectrum has a simple pole at g = gi. [33] In the finite chain 
with N = 10, 000, the singularity is replaced by a very sharp maximum. Introduction of even 
slight disorder tends to further broaden and randomize this peak. Obviously, this broaden- 
ing, as well as the randomization of the spectrum in the vicinity of g = gi, result in spatial 
decay of the ordinary SP excitated by a near-field tip. However, the peak corresponding 
to the extraordinary SP is almost unaffected by disorder. Consequently, the extraordinary 
SP does not experience localization and related spatial decay when off-diagonal disorder is 
introduced into the system. 

The physical reason why the extraordinary SP is not affected by disorder is quite straight- 
forward. This SP is mediated by electromagnetic waves in the far (radiation) zone that arrive 




(15) 
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at a given nanosphere from all other nanospheres. The synchronism condition (that all these 
secondary waves arrive in phase) is not affected by disorder as long as the displacement am- 
plitudes are small compared to the wavelength. In the numerical examples of this section, 
wavelength is ten times larger than the inter-particle spacing, which is, in turn, much larger 
than the displacement amplitudes. Additionally, effects of disorder are expected to be av- 
eraged if the amplitudes are mathematically independent because the scattered field at 
a given nanosphere is a sum of very large number of secondary waves. In contrast, the 
ordinary SP is mediated by near-field interactions which are very sensitive to even slight 
displacements of nanospheres. Besides, the scattered field at a given nanosphere (when the 
ordinary SP propagates in the chain) is a rapidly converging sum of secondary fields, so that 
only a few terms in this sum are important and no effective averaging takes place. 

B. Diagonal Disorder 

Diagonal disorder is disorder in the properties of nanospheres. There are several possi- 
bilities for introducing such disorder. 

First, the spheres can be polydispcrsc, i.e., have different radiuses a„. This is, however, 
not a truly diagonal disorder. Indeed, polarizability of n-th nanosphere can be written in 
this case as /„(q;) where = (a„/ (a))^ and (a) is the average radius (note that the radiative 
correction to l/ctn can be included into the dipole sum S{k, q), so that this analysis remains 
valid even when this correction is important). The factors /„ are all positive-definite which 
allows one to introduce a simple transformation of Eq. (1) which removes the diagonal 
disorder [34]. Then the disorder becomes effectively off-diagonal. More importantly, one 
can retain a well-defined spectral parameter of the theory, l/{a). We, therefore, do not 
consider polydispersity in this section. 

The second possibility is variation of the absorptive parameter 7. This effect that can be 
practically important. Yet, we are interested in propagation in the absence of Ohmic losses 
and, therefore, set 7„ = 0. 

The third, and the most fundamental, reason for diagonal disorder is variation of the 
Frohlich frequency of nanospheres. Namely, we take the Frohlich frequency of the n-th 
particle to be u-pn = (^f)(1 +'Cn)) where ^„ are statistically independent random variables 
evenly distributed in the interval [— A, ^4] and uj/{uj-p) — 1- The factors /„ in this case are 
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no longer positive-definite and the transformation of Ref. 34 can not be applied. The most 
profound consequence of introducing the diagonal disorder is that the spectral parameter 
such as l/a is no longer well defined. As long as the disorder amplitude is relatively small, 
one can view 1/(q;) as an approximate spectral parameter. However, as the amplitude of 
disorder increases, this approach becomes invalid. We have seen in section IV that variation 
of the ratio uj/uj-p in the interval 0.988 < uj/uj-p < 1.007 can result in dramatic changes in the 
way an SP excitation propagates along the chain. For a;/a;F outside of this interval, ordinary 
SP could not be effectively excited. However, in section IV, variation of lo/lof apphed to 
all nanospheres simultaneously. We now introduce random uncorrelated variation of this 
ratio for each individual nanosphcrc. The effects of such disorder are difficult to predict 
theoretically. We can, however, expect that these effects become dramatic for A > 0.01 
since, in this case, the ratio uj/uj-pn can be outside of the interval [0.988, 1.007]. 

We first adduce in Fig. 10 the same dependencies as shown in Fig. 5 but for different levels 
of diagonal disorder. The results appear to be, qualitatively, quite similar, although in the 
case A — 0.016, there is no visible trend lor x/h > 3000. We then investigate whether the 
extraordinary SP remains insensitive to diagonal disorder. Data analogous to those shown in 
Fig. 8, but for diagonal disorder, arc presented in Fig. 11. It can be seen that the influence of 
diagonal disorder on the extraordinary SP is stronger than that of the off-diagonal disorder. 
When the disorder amphtude A exceeds 0.01, the influence becomes quite dramatic. But 
paradoxically, A — 0.016, the average decay rate is much slower than for A — 0.008, 
although the amplitude of fluctuations is much larger. 

To see why this happens, we look at the extinction spectra ae{q). These are plotted 
in Fig. 12. The fundamental difference between the off-diagonal and diagonal disorder is 
clearly revealed by comparing this figure to Fig. 9. Namely, the effect of diagonal disorder 
is not only to broaden and randomize the peak aX q — q^, but also to shift it towards the 
peak corresponding to the extraordinary SP. At sufficiently large levels of disorder [A > 
0.016), the separate peaks disappear and a broad structure emerges. At this point, ordinary 
and extraordinary SP can no longer be distinguished. Correspondingly, the ordinary and 
extraordinary SP are effectively mixed in the A = 0.016 curve in Fig. 11, while fo smaller 
amplitudes of the diagonal disorder, the extraordinary SP is excited predominatUy. 
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VI. SUMMARY 



We have considered surface plasmon (SP) propagation in a linear chain of metal nanopar- 
ticles. Computer simulations reveal the existence of two types of plasmons: ordinary (qua- 
sistatic) and extraordinary (non-quasistatic) SPs. 

The ordinary SP is characterized by short-range interaction of nanospheres in a chain. 
The retardation effects are inessential for its existence and properties. The ordinary SP 
behaves as a quasistatic excitation. The ordinary SP can not radiate into the far zone in 
perfectly periodic chains because its wave number is larger than the wavenumber k — uj/coi 
free electromagnetic waves. However, it can experience decay due to absorptive dissipation 
in the material. 

The second, extraordinary, SP propagates due to long-range (radiation zone) interaction 
in a chain. Its excitation is possible due to the existence of the non-Lorentzian optical 
resonance in the chain introduced in Ref. 16. The extraordinary SP may experience some 
radiative loss but is much less afected by absorptive dissipation and disorder. As a result, it 
can propagate to much larger distances along the chain. The extraordinary SP can be used 
to guide energy or information in all-optical integrated photonic systems. 

We have also considered the effects of disorder and localization of the ordinary and 
extraordinary SPs. Results of numerical simulations suggest that even small disorder in the 
position or properties of nanoparticles results in localization of the ordinary SP. However, the 
extraordinary SP appears to remain delocalized for all types and levels of disorder considered 
in the paper. 

The authors can be reached at vmarkel@mail.med.upenn.edu (VAM) and 
asarychev@ethertronics.com (AKS). 
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FIG. 10: Same as in Fig. 5, but for different levels A of diagonal disorder. 
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FIG. 11: Same as in Fig. 10 but for uj/oj-p = 0.984 and for different levels A of off-diagonal disorder, 
as indicated. Note that the curves for A = and A = 0.01 are indistinguishable. 
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FIG. 12: Specific extinction Ce as a function of lateral wave number of incident wave, g, for different 
levels of off-diagonal disorder, A. 
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